Two-dimensional solitons and vortices in media with incommensurate linear and 

nonlinear lattice potentials 



Jianhua Zeng 1 ' 2 and Boris A. Malomed 1 

1 Department of Physical Electronics, School of Electrical Engineering, 
Faculty of Engineering, Tel Aviv University, Tel Aviv 69978, Israel 

2 State Key Laboratory of Low Dimensional Quantum Physics, 
Department of Physics, Tsinghua University, Beijing 100084, China 

We construct families of ordinary and gap solitons (GSs), including solitary vortices, in the two- 
dimensional (2D) system based on the nonlinear-Schrodinger/Gross-Pitaevskii equation with the 2D 
or quasi- ID (Q1D) periodic linear potential, combined with the periodic modulation of the cubic 
nonlinearity (also in the 2D or Q1D form), which is, generally, incommensurate with the linear 
potential, thus forming a "nonlinear quasicrystal" . Stable vortices are built as complexes of four 
peaks with the separation between them equal to the double period of the linear potential. The 
system may be realized in photonic crystals or Bose-Einstein condensates (BECs). The variational 
approximation (VA) is applied to ordinary solitons (residing in the semi-infinite gap), and numerical 
methods are used to construct solitons of all the types. Stability regions are identified for soliton 
families in all the versions of the model. 

PACS numbers: 05.45.Yv, 03.75.Lm, 42.65.Tg 

I. INTRODUCTION 

A versatile technique for the control of guided photonic and matter waves is based on the use of periodic (lattice) 
potentials. In the case of Bose-Einstein condensate (BEC) the potential may be induced by optical lattices, and in 
photonics — by transverse structures in the form of photonic crystals |l| . In nonlinear media, the lattice potentials help 
to create and stabilize various types of solitons, both ordinary ones (found in the semi-infinite gap of the underlying 
spectral structure [2]) and gap solitons (GSs) which exist in the presence of the self-attractive and repulsive 
nonlinearity, respectively. Along with the fundamental solitons, the periodic potentials may support stable solitary 
vortices, built as multi-peak patterns with the phase distribution carrying the topological charge, in the case of the 
self- attraction [2[ and self-repulsion |4j alike. 

Results accumulated in theoretical and experimental studies of one- and two-dimensional (ID) and (2D) solitons 
(including 2D vortices) supported by linear-lattice (LL) potentials were reviewed, respectively, in Refs. [Ej and [fj, 
see also a more recent review A related topic is the study of discrete solitons in optics, which correspond to the 
limit of very deep lattice potentials 0|. 

The photonic-crystal structures induce, simultaneously with the LL potentials, an effective nonlinear potential (alias 
pseudopotential [9]), in the form of the concomitant spatially periodic modulation of the local nonlinearity coefficient. 
In BEC, nonlinear lattices (NLs) may be induced by external fields which affect the local nonlinearity via the Feshbach 
resonance. The studies of solitons in NLs, as well as in LL-NL combinations, have been recently reviewed in Re f. flQj . 
A direct experimental observation of NL-supported optical solitons was reported at a surface between lattices . 

A natural generalization of the setting combining the LL and NL is the one with different or incommensurate 
periodicities of the two lattices. In the framework of the ID setting, both ordinary solitons and GSs, supported by 
such a combination of competing linear and nonlinear potentials, were studied in Ref. [12| by means of numerical 
methods and analytical approximations. Noteworthy results were obtained, in particular, for existence borders of the 
solitons as functions of the LL-NL incommensurability, and for the empirical "anti-Vakhitov-Kolokolov" (anti-VK) 
stability criterion for GSs, which is written in terms of the dependence of the chemical potential, /i, on norm, A/", of the 
soliton: dfi/dN > (the VK criterion for ordinary solitons supported by the self- attraction in the semi- infinite gap is 
dfi/dN < [HI). The objective of the present work is to produce results for 2D solitons supported by incommensurate 
LL-NL combinations. We develop a variational approximation (VA) which is applied, along with numerical methods, 
to ordinary solitons, while GSs in the first finite bandgap are studied solely in a numerical form, as the analytical 
approach would be too cumbersome in that case. Numerical methods are also used to construct vortex solitons, in 
the semi-infinite and finite gaps alike. 

The settings considered here include both full 2D lattices potentials and quasi- ID (Q1D) ones, which depend on 
the single coordinate (the LL potential of the Q1D type is sufficient for the stabilization of ordinary 2D solitons, in 
diverse realizations of 2D media with the self- attractive nonlinearity [14(). In fact, the combination of the periodic but 
mutually incommensurate linear and nonlinear lattices makes the medium effectively tantamount to a quasicrystal 
for nonlinear excitations. Fundamental solitons and solitary vortices in linear quasiperiodic potentials were studied 



FIG. 1: The band structure of the linearized equation (pQ) in the first Brillouin zone, for depth e = 7.4 of the linear lattice; \i 
is the chemical potential of the Bloch waves. 



theoretically [15|, and 2D pho tonic quasicrystals have been recently created experimentally [lfj. It is also relevant 
to mention a recent work [17], which was dealing with 2D solitons in a model combining crossed Q1D linear and 
nonlinear periodic potentials. 

The rest of the paper is organized as follows. The model is introduced in Section II, ordinary and gap solitons are 
considered, severally, in Sections III and IV (each section reports the results for fundamental and vortex solitons), 
and the work is concluded by Section V. 



The 2D system combining the periodic LL potential and NL pseudopotential may be written in the form of the 
scaled Gross-Pitaevskii (or nonlinear Schrodinger) equation for the BEC wave function (or the local amplitude of the 
electromagnetic wave guided by the photonic crystal), ip (x,y,t) [1(1 Il2j: 



where t is time (or the propagation distance in the photonic-crystal waveguide), Laplacian V 2 = <9 2 + <9 2 acts on 
transverse coordinates x and the scale in the (x,y) plane is set by fixing the LL period to be 1, the period of the 
NL is 2/q (that is, q is the incommensurability index), and the NL strength is normalized to be g = ±1. The center 
of the soliton will be placed at point x = y = 0, hence g = +1 and —1 correspond, severally, to the dominating 
self-attraction or self-repulsion, which support ordinary solitons in the semi- infinite gap, or GSs in finite bandgap(s), 
respectively. 

The remaining parameter in Eq. (pQ) is the normalized LL strength, e. Generic results for the case when the 
system's spectrum contains a single finite bandgap are reported here for a fixed strength, e = 7.4. The respective 
band structure in the first Brillouin zone [1], found from the linearized version of Eq. (pQ), is displayed in Fig. [TJ 

At q = 0, Eq. (pQ) amounts to the usual model with the LL potential and spatially uniform nonlinearity. The 
LL and NL are commensurate at q = 2, and the subharmonic commensur ability, with the ratio of the LL and NL 
periods 1:2, occurs at q = 1. The full incommensurability (overall quasi-periodicity in the system) corresponds to 
irrational values of q, but, in practical terms, the incommensurability may be emulated by q = 1.5, with the period 
ratio 3 : 4. The Q1D versions of the model correspond to dropping terms cos (27n/) and/or cos (irqy) in the linear 
and/or nonlinear potentials. 

Stationary solutions to Eq. (pQ) with chemical potential \i (or propagation constant — /i, in terms of the guided 
optical waves) are looked for as ip(x,y,t) = <j>(x, y) exp(— ifit) : where function (j)(x,y) satisfies equation 



II. THE MODEL 



iipt = -(l/2)V 2 ^-^[cos(27rx)+cos(27r?/) ] ^ 
— g [cos(irqx) + cos(irqy)] \ip\ 2 ^. 



(1) 



/i0 



(l/2)V 2 (/> - e [cos(2ttx) + cos(27n/) ] <j> 
g [cos(7rgx) + cos(iiqy)\ |(/>| 2 (/>. 



(2) 



The VA will be based on the Lagrangian of Eq. (j2j), which is 
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FIG. 2: (Color online) The chemical potential vs. the norm for soliton families in the semi-infinite gap of the model with the 
full 2D linear and nonlinear potentials, as produced by the VA (dashed curves), and as found from numerical solutions of Eq. 
([2]) (chains of symbols). 




FIG. 3: (Color online) A typical example of a stable fundamental soliton found in the semi-infinite gap of the model with the 
full 2D linear and nonlinear potentials, for # = +1, #=1.5, and \i — —5.24, N = 1.22. The right panel shows contour plots of 
the stationary real wave function, <j> (x, y), for the soliton. 



III. LOCALIZED MODES IN THE SEMI-INFINITE GAP 



A. Fundamental solitons 



We start the analysis of the ordinary solitons, which are expected to exist in the semi- infinite gap at g = +1 in 
Eq. (pQ), using the VA based on the Gaussian ansatz, cj) (x,y) = Aexp [— (x 2 + y 2 ) / (2W 2 )] , with the corresponding 

norm N = J J<p 2 (x, y) dxdy = it (AW) 2 @. The substitution of the ansatz into the Lagrangian yields the following 
expression, written in terms of N and width W: 



T N 
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and the respective variational equations, dL e g/dW = dL e g/dN = 0: 
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Dependences fi(N) for the soliton families, produced by a numerical solution of Eqs. ((4]) at different values of 
incommensurability index q, are displayed in Fig. [2j along with their counterparts, produced by numerical solutions 
of stationary equation (j2j) and verified in direct simulations of Eq. (PQ) . 

Shapes of the stable solitons generated in the semi-infinite gap by the full 2D model, and by its version with the Q1D 
linear potential, are displayed in Figs. [3]and|4l Naturally, these shapes are, respectively, quasi- isotropic and strongly 
elongated, resembling those reported previously in other 2D models stabilized by the LL potentials (2I. Irl HoL fl7l| . 

It is seen from Fig. [2] that the VA describes the ordinary solitons with a reasonable accuracy, except near the 
edge of the semi-infinite gap, where the Gaussian ansatz is irrelevant, as the soliton becomes very wide and develops 
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FIG. 4: (Color online) The same as in Fig. [3j but in the case of the Q1D linear potential (while the nonlinear potential remains 
two-dimensional), for g = +1, q = 1, and \i — —2.91, N = 1.68. The shape of the solitons in the version of the model where 
the nonlinear potential is also made quasi-one-dimensional is quite similar to the present one. 
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FIG. 5: (a) Stability borders for the family of the ordinary solitons (in the semi-infinite gap) in the model with the 2D linear 
potential and 2D or Q1D nonlinear one (solid and dashed lines, respectively), (b) The same, in the case of the Q1D linear 
potential. The solitons are stable beneath the respective borders. 



a complex shape. Further, simulations of the evolution of perturbed solitons demonstrate that the stability of the 
solitons exactly obeys the VK criterion, d/a/dN < [strictly speaking, if dependence fi(N) is taken in the numerically 
found form]. These features of the families of ordinary-soliton solutions are similar to those found before in the ID 
variant of the model [12j. 

The stability of the 2D solitons in the semi-infinite gap, for all the four realizations of the model (2D or Q1D linear 
and nonlinear potentials) is summarized in Fig. [5] When solitons are unstable, they suffer decay into radiation waves 
(rather than rearranging themselves into stable solitons). It is worthy to note that the replacement of the full 2D 
NL by its Q1D counterpart leads to an expansion of the stability areas, which is explained by the fact that, in the 
case of the Q1D NL, the ordinary solitons should make the effort to "dodge" the destabilizing locally self-repulsive 
nonlinearity only in one direction, rather than in two. 

In addition to the above analysis, we tried to test the mobility of the solitons by simulating their evolution after 
sudden application of a "kick", i.e., multiplication of the wave function of a stable quiescent soliton by the phase-tilt 
factor, exp (i (k x x + k y y)), with vectorial kick parameter k. Except for the obvious case when both the LL and NL 
have the collinear Q1D structure, and the kick is applied in the unconfined direction, mobile solitons were not found, 
even if either the LL or NL potential was of the Q1D type. Instead of setting the soliton in motion, a sufficiently strong 
kick tends to destroy it. 



B. Solitary vortices 



Stable solitary vortices with topological charge 1 were found as "hollow" four-peak complexes, with the separation 
between the peaks equal to the double period of the LL potential (Ax = Ay = 2), and a nearly empty cell at the 
center of the pattern, see an example in Fig. [6] It is relatively easy to find stable vortices of this type, due to the 
weak interaction between the peaks. More densely packed vortex patterns can also be constructed, but we could not 
find stability regions for them. It is known from the studies of other models too that the vortices with inner "voids" 
are more likely to be stable 0, [To| . 



FIG. 6: (Color online) An example of a stable vortex with topological charge 1, which is supported, in the semi- infinite gap, 
by the combination of the full 2D linear and nonlinear structures with q = 1 and g — +1. The left, middle, and right panels 
display, respectively, the absolute value of the field, \(f>(x,y)\, as a function of the coordinates, contour plots of \(j)(x,y)\, and 
the distribution of the phase, which carries the vorticity. Parameters of the vortex are fi — —5.1, N — 4.4. 




FIG. 7: (Color online) The same as in Fig. 2, but for families of solitary vortices in the semi-infinite gap, supported by the 
combination of the full 2D linear potential with the 2D (a) or Q1D (b) nonlinear lattice. 



Actually, the vortices are stable only for values of the incommensurability index close to q = 0, 1, and 2 — namely, 
within intervals of half- width Aq ~ 0.1 around these values. The latter observation may be explained by the fact 
that, at such values of both the linear and nonlinear potentials have minima at or close to sites where the the power 
(density) peaks are located. 

Furthermore, the vortices of the type shown in Fig. [6] are found to be stable (in the semi-infinite gap) in the case 
when the LL potential is fully two-dimensional, while the NL may be of either 2D or Q1D type. The corresponding 
families of the vortex modes are represented by fi(N) curves in Fig. [71 Detailed analysis demonstrates that the 
stability of these families exactly follows the VK criterion, i.e., stable are portions of the families with dji/dN < 0. 



IV. GAP SOLITONS AND VORTICES 



A. Fundamental solitons 



Numerically generated fundamental GSs feature, as usual, more complex shapes than the ordinary solitons, see an 
example of a stable GS, found sufficiently far from edges of the finite bandgap, in Fig. [8] (in the model combining 
the 2D LL and Q1D NL, the shapes of the GSs are quite similar). The GS families are characterized by /jl(N) curves 
which are shown in Fig. [9]^a). Unlike the ordinary solitons (cf. Fig. [2]), the GSs always feature dfi/dN > 0, thus 



FIG. 8: (Color online) A typical example of a stable gap soliton in the model with the full 2D NL, for g = —l,q = 1.5, and 
H — 3.4, N — 1.4. 




FIG. 9: (a) Curves /jl(N) for families of gap solitons in the model with the full 2D nonlinear potential. Black squares denote 
unstable solitons close to the bottom edge of the bandgap. (b) Stability borders for the entire set of the GS families in the 
model with the 2D and Q1D nonlinear potentials (solid and dashed lines, respectively). The solitons are stable between the 
corresponding stability borders. Both the top and bottom borders are located near the respective edges of the finite bandgap. 



complying with the "anti-VK" criterion proposed in Ref. [12 ]. 



Numerical tests demonstrate that the GSs tend to be stable sufficiently deep inside the finite bandgap, and unstable 
near its edges (unstable GSs suffer decay into radiation). The numerically found stability borders for the entire set 
of GSs in the two versions of the present model, with the full 2D NL and its Q1D counterpart, are presented in Fig. 
Efb). It is observed that, on the contrary to the solitons in the semi-infinite gap, the stability region for GSs tends 
to be essentially narrower for the NL of the Q1D type, in comparison with the full 2D NL. The latter feature seems 
natural, as, unlike the case of ordinary solitons, the NL may provide for an additional support to the GSs. 

With the increase of the GS stability areas clearly tend to shrink to nil, which actually happens in Fig. [9f b) for 
the variant of the model with the Q1D NL (we expect the same ought to happen for the full 2D NL, but numerical 
problems impede extending the stability diagram to still larger values of q). This trend is explained by the fact that, 
at large the rapidly oscillating NL field tends to average itself to zero, hence the broad (see Fig. [8|) GS ceases to feel 
the action of the nonlinearity. The ordinary solitons in the semi-infinite gap do not demonstrate such a trend (cf. Fig. 
Ej), as, following the increase of q, these solitons are able to compress themselves inside a single cell of the structure, 
remaining centered around a region with the self- attractive nonlinearity. Finally, numerical tests demonstrate that, 
as well as it was concluded above for the ordinary solitons, the GSs are not mobile objects (not shown here in detail). 



B. Solitary vortices 



Stable vortex solitons that can be found in the finite bandgap feature the same structure which supports the stable 
vortices in the semi-infinite gap (cf. Fig. [6]), i.e., they are built of four peaks separated by the distance equal to the 
double LL period, the vorticity being carried by the corresponding phase distribution, see an example in Fig. [TUJ The 
difference from the situation in the semi-infinite gap is that the solitary vortices in the finite bandgap may be stable 
only when both LL and NL have the full 2D structure (i.e., the vortices are unstable if the NL is of the Q1D type). 

As well as their counterparts in the semi-infinite gap, the solitary vortices in the finite bandgap are found to be stable 
only for the values of incommensurability index close to q = 0, 1, and 2. Families of these vortices are represented in 
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FIG. 10: (Color online) An example of the stable vortex soliton found in the finite bandgap, for q— 1, g — —1 and \i — 3.5, N = 
5.36. The meaning of the three panels is the same as in Fig. [6j 




N 



FIG. 11: The dependence of the chemical potential on the norm for families of solitary vortices in the finite bandgap. Stable 
portions of the vortex families are located inside the marked stripe. 

Fig. E] by the corresponding fi(N) curves. As well as the fundamental gap solitons, the stable vortices in the finite 
gap obey the "anti-VK criterion, dfi/dN > 0. 

V. CONCLUSIONS 

We have introduced the model of 2D nonlinear photonic crystals and BEC based on the interplay of linear and 
nonlinear lattices with different (generally, incommensurate) periods, which may be considered as a "nonlinear qua- 
sicrystal". Both fully 2D periodic potentials and their Q1D reductions were considered. For the ordinary solitons in 
the semi-infinite gap, the VA (variational approximation) was developed. In the general case, the solitons and solitary 
vortices were explored by means of numerical methods. The stability regions have been identified for the entire sets 
of the ordinary solitons and GSs (gap solitons). Stable families of vortex solitons of both types have been found too. 

A relevant direction for the development of the analysis reported above may be search for stable higher-order 
vortices. On the other hand, it may also be interesting to extend the analysis to a broader parameter region, which 
may give rise to higher-order bandgaps, in addition to the single finite bandgap existing in the situation considered 
in this work, and construct solitons and solitary vortices in the additional gaps. 
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